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ABSTRACT 

The modal element method has been employed to determine the 
scattered field from a plane acoustic wave impinging on a two 
dimensional body. In the modal element method, the scattering body 
is represented by finite elements, which are coupled to an 
eigenfunction expansion representing the acoustic pressure in the 
infinite computational domain surrounding the body. The present 
paper extends the previous work by developing the algorithm 
necessary to calculate the acoustic scattering cross section by the 
modal element method. The scattering cross section is the acoustical 
equivalent to the radar cross section (RCS) in electromagnetic 
theory. Since the scattering cross section is evaluated at infinite 
distance from the body, an asymptotic approximation is used in 
conjunction with the standard modal element method. For valida- 
tion, the scattering cross section of the rigid circular cylinder is 
computed for the frequency range 0.1 < ka < 100. Results show 
excellent agreement with the analytic solution. 


INTRODUCTION 

The modal element method, which couples finite element algo- 
rithms to eigenfunction expansions, has been employed in calculat- 
ing the scattered field from an acoustical plane wave impinging on 
a two dimensional body. The primary reasons for employing the 
modal element method are (1) to accurately describe the radiation 
boundary condition at the computational boundary and (2) to reduce 
the size of the numerical grid. In fact, for hard scatterers, the modal 
element method can effectively reduce a two dimensional scattering 
problem to a one dimensional problem by employing a single line of 
elements circumscribing the scattering body. 

The modal element method has been given various titles, such as 
the unimoment method and the transfinite element method. In 


electromagnetics, the method has been applied to scattering from 
dielectric cylinders (Chang & Mei (1976), Lee & Cendes (1987), 
Baumeister & Kreider ( 1 992)) and propagation in ducts (B aumeister 
(1991 )). In acoustics, the method has also been applied to scattering 
from cylinders (Baumeister & Kreider (1993)) and propagation in 
ducts (Astley & Eversman (1981)). In all of these applications, an 
eigenfunction expansion is used to represent the acoustic pressure 
field in the far field. An asymptotic approximation of this expansion 
presents a simple means of determining the scattering cross section. 

This paper presents the numerical algorithm for evaluating the 
acoustic scattering cross section by the modal element method. The 
scattering cross section is the acoustical equivalent to the radar cross 
section (RCS) in electromagnetic theory. Since the scattering cross 
section is evaluated at infinite distance from the body, asymptotic 
approximations are used in conjunction with the standard modal 
element method. For validation, the method is applied to scattering 
from rigid circular cylinders, for which the analytic solution is 
known. 


NOMENCLATURE 

A* modal amplitude of wave moving radially outwards 
a dimensionless circular cylinder radius 

H (1) Hankel function of the first kind 

m 

k wave number 

m mode number 

M number of modal coefficients used in eigenfunction expansion 

cocf 
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p dimensionless perturbation acoustic pressure 

r dimensionless radial coordinate 

e dimensionless property constant 

0 angle between radius vector and positive x axis 

p dimensionless property constant 

o c acoustic scattering cross section 

co dimensionless frequency 

Superscript 

* complex conjugate 


In the finite element domain, an approximate solution for the total 
(incident + scattered) pressure at the element nodes is calculated by the 
Galerkin method. Linear triangular elements are used and the subdomain 
interface is approximated by piecewise linear segments. In the homo- 
geneous domain, which extends to infinity, the pressure is represented 
by an eigenfunction expansion. The modal element method couples 
the two solution forms by imposing continuity on the pressure and 
velocity at the interface between the two subdomains. This coupling 
results in a single matrix equation in which the eigenfunction coeffi- 
cients and the pressure at the finite element nodes are calculated 
simultaneously, yielding a global representation of the acoustic pres- 
sure field. 


GOVERNING EQUATION AND BOUNDARY CONDITIONS 

The scalar form of the acoustic wave equation can be written as 
(Baumeister & Kreider (1993)) 


METHOD OF ANALYSIS 

The present study is concerned with computing the scattering cross 
section of a two dimensional axi symmetrical rigid body due to an 
impinging plane wave traveling in the +x direction. In order to 
determine the acoustic scattering cross section, the scattered field must 
be computed. To accomplish this, the spatial domain is divided into 
two subdomains, the homogeneous domain and the finite element 
domain, as shown in Fig. 1 . The grid system in Fig. 1 is designed to 
allow various structures to be imbedded into the grid, as shown by 
inserts A and B. For the special case of rigid bodies, the grid system can 
be reduced to a single ring of elements, which effectively reduces the 
two dimensional problem to a one dimensional problem. 
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Figure 1 . — Finite element grid system. 



+ co 2 jxp = 0 


( 1 ) 


The harmonic time dependence e' iux has been factored out. e and p are 
material property constants and to is the dimensionless frequency. The 
wave number is 

k = -J co 2 pe (2) 


In the homogeneous region surrounding the body, e= 1 and p=l , while 
inside the body, e and p may assume other appropriate values depend- 
ing on the material. Baumeister and Dahl (1989) present explicit 
expressions for e and p as a function of medium density, porosity, 
viscous loss coefficient and a heat transfer parameter. 

At the interface between the finite element region and the analytic 
region, continuity is imposed on the pressure and velocity. The 
radiation boundary condition at infinity is automatically satisfied by 
the eigenfunction expansion introduced in the next section. 


ANALYTIC SOLUTION 

In the homogeneous domain, an exact expansion for the pressure can 
be derived (Morse and Ingard, 1 968, p. 401 ) from eq. ( 1 ) by separation 
of variables: 


_ _ ikx . 

p = p+p =e + 


COCI 

X A m H m (kr)COS(m0) - 


(3) 


m=0 


p j is the incident plane wave; p s is the scattered wave. This expansion 
is valid for symmetric bodies. The modal coefficients A m + are un- 
knowns to be determined. Formulas to estimate the number M f of 
modes needed for convergence can be found in Baumeister & Kreider 
(1993). 
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FINITE ELEMENT SOLUTION 


The solution to eq. (1) will be obtained by linear finite element 
theory. In Fig. 1, the finite element region is divided into triangular 
elements defined by comer nodal points. It is assumed that all material 
properties are constant in each element. In the conventional weighted 
residual approach, the unknown pressure field is described in terms of 
all the unknown nodal values of the pressure. 

The finite element aspects of converting eq. (1) and (3) and bound- 
ary conditions into an appropriate set of global difference equations 
can be found in Baumeister & Kreider (1993). The resulting set of 
global difference equations is solved for the pressure at the nodes and 
the modal coefficients A m + . With the modal coefficients determined, a 
relatively simple algorithm has been developed to determine the 
scattered cross section and is described in the next section. 


o c (0) = lim 27cr| 

r— ><» 


I o if kr-^-n V“coet * _i m — 

J— 4 ' Y, A m e 2 cos(m0) 
inkr m 

m=U 


M f -1 . * 

coet -im— 



' M coef-> 


K 

im— 


^ -ifkr-— Jt] __ 11W 

1 e ^ 4 ^ 2 cos(m0) 

m=0 




J 


^T 1 _ im n 


YM «rf“l im— 


''-m — till — v— ^ * tan 

2 cos(m0) | ^A*e 2 cos(m0) 

^ m=0 


m=0 


(7) 


BISTATIC ACOUSTIC SCATTERING CROSS SECTION 

The bistatic acoustic scattering cross section per unit length for a two 
dimensional body can be defined as 

|p s | 2 (p s «p s *) 

O c (0)= lim 2m - — y = lim 27irj— — r (4) 

r "»“ |p'| T ~*°° (P *P ) 

The standard superscript c implies a two dimensional problem, where 
the cross section is often called the “echoing width” (Ruck, vol. 1, 
pp. 23). In three dimensions, the scattered field decays as 1/r 2 while in 
two dimensions, considered herein, the scattered field decays as 1/r. 
Thus, the 27tr is employed in eq. (4) to negate the 2D scattered field 
decay and to yield a finite measure of the reflected power at infinitey. 
The cross section defined in eq. (4) could be used in range equations 
specifically defined for 2D bodies. 

For a unit incident plane wave, eq. (4) simplifies to 

cr c (0)= lim 2jtr(p s »p s *) (5) 

r->oo ' ' 

Using conventional numerical grids, it is often necessary to extrapo- 
late the nodal values of the field on the outer boundary to infinity. For 
example, Noack and Anderson (1992) include a boundary element 
approximation to determine the scattering cross section at infinity. 
Fortunately, no such additional work is required with the modal 
element method. Since the scattering cross section is evaluated at 
infinity, the asymptotic expansion forthe Hankel functions (Abramowitz 
& Stegun, 9.2.3) 



can be inserted into eq. (3), which is then substituted into eq. (4) 
to yield 


Once the A + have been calculated, the scattering cross section and the 

m 

pressure field can be easily determined. 

RESULTS AND COMPARISONS 

Consider the following scattering problem: a unit plane wave, 
incident from the left, strikes a rigid circular cylinder of dimensionless 
radius a = 1 oriented with its axis normal to the propagation direction. 
The rigid body is simulated numerically through an impedance mis- 
match induced by setting £ = 1 - 1 0 19 i and |jl = 1 for each internal finite 
element and £ = 1 and \i = 1 in the homogeneous domain surrounding 
the body. This feature allows flexibility in the numerical implementa- 
tion of the method. Rigid bodies with absorber coatings may be studied 
with only slight modification to the computer code by adding several 
finite element rings with the appropriate values of £ and \i. 

This problem is solved for frequencies ranging from ka = 0.1 to 
ka = 100 in order to validate the method. This problem was chosen 
because it has an exact solution, given by eq. (3) with modal coeffi- 
cients (Bowman et al., eq. (2.38)) 



and because it has been examined using other methods (Ling, 1986, 
Fig. 15). 

In Fig. 2, the computed cross sections (square symbols) are com- 
pared to the corresponding exact solutions (solid line). The results 
agree very well with the exact solutions over a wide frequency range. 
Results for ka = 50 and ka = 100 show similar agreement but are 
omitted for brevity. In addition, the results agree well with the low 
frequency calculations (k = 3.2) of Ling (1986, Fig. 15). 

Additional validation of the algorithm is obtained by examining 
high frequency (ka > 20) asymptotic approximations. In the 
backscattered direction (0 = 180), the optical approximation is valid 
(Ruck, et al. (1970), eq. (4.1) to (38)) and the backscattering cross 
section becomes 

a c = rca = 7t (a - 1) (9) 
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Figure 2. — Bistatic scattering cross section. 


In the forward scattered direction (0 = 0), commonly referred to as the directions. Even at ka = 1 00, the relative error in the forward direction 

diffracted field, the optical approximation yields (Ruck, et al. (1970), is only about 3 percent. 

eq. (4. 1) to (41 )) the forward cross section Some error at high frequencies can be attributed to the approxima- 

tion involved with the asymptotic expansion given by eq. (7). How- 

a c - 4ka 2 - 4k (a - 1 ) (10) ever ’ ma J or source °f error l ^ e calculated values of the mode 

coefficients A m + ineq. (3). For example, in the k= 100 calculation, over 

130 modal coefficients need to be determined. Consequently, some 
Figure 3 shows that the asymptotic limits are attained as frequency significant errors m involved in the caJcu i at i 0 „ 0 f the higher order 
increases in both the forward and backscattered directions. Table I mo d es 
contains the numerical values obtained for various frequencies in these 
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Forward and reverse cross sections 



TABLE I. — COMPARISON OF EXACT AND 
MODAL ELEMENT (ME) BISTATIC 
CROSS SECTIONS 


k 

Forward (0 = 0) 

Backward (0 = 180) 


Exact 

ME 

Exact 

ME 

0.1 

0.00261 

0.00263 

0.0022 

0.00219 

1.0 

1.643 

1.659 

3.427 

3.437 

3.2 

8.459 

8.470 

3.412 

3.413 

5.0 

14.77 

15.22 

2.795 

2.778 

25 

90.31 

91.00 

3.151 

3.153 

50 

187.6 

189.3 

3.140 

3.140 

100 

378.7 

369.7 

3.106 

2.957 


CONCLUDING REMARKS 

In the modal element method, the scattering body is imbedded in a 
finite element grid to calculate the acoustic pressure field. In the far 
field, the pressure is represented analytically by a Hankel function 
expansion. The asymptotic approximation of the Hankel functions are 
used to calculate the scattering cross section. The method is applicable 
to problems involving very high or low frequency scattering from 
bodies. In the validation examples, the numerical results are in good 
agreement with the corresponding exact solutions and asymptotic 
analytical limits. 
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